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I. Introduction 

This FORTRAN program produces Calcomp plots of selected, 
portions of the celestial sphere. One example is shown in 

Figure 1. The brightest stars in the regions are selected 

/ 

out of two magnetic tapes which are rewritten versions of 
the SAO Star Catalog tapes. The stars are plotted in their 
correct positions relative to the center for a flat plate 
photograph, and are indexed in a table with further informa- 
tion on both plot and line printer. Positions of x-ray sources 
can be input and will be plotted and listed, and any nearby 
stars, irrespective of magnitude, can be found. Also, a 
variable mesh grid of constant a and S lines and a title can 
be specified. Stars may be selected by spectral type. 

II. Directions for Using the Program 

To use -the program, preface the deck*with the standard 
JCL cards, and follow it with the following: 


//G.FT01F001 DD UNIT® (TAPE9 ,, DEFER) , LABEL® ( ,NL) , + 

// V0LUME=SER= XR4 026 ,DISP= (OLD, PASS) , + 

// DCB® ( RECFM= V , LRECL=7 188, BLKS I ZE= 7 2 0 0 ) 

//G . FTO 2F0 01 DD UNIT- (TAPE 9 , , DEFER) , LABEL= ( ,NL) , + 

// VOLUME=SER= XR4Q27 ,DISP= (OLD, PASS) , + 

// DCB® ( RECFM=V , LRECL® 89 88, BLKS I ZE=9 0 0 0 ) 

//G . PLOTTAPE DD UNIT® (TAPE 7 , , DEFER) , LABEL® ( , BLP ) , * 

// V0LUME=SER= XR2 084 

//G. FT05F002 DD DUMMY 
//G. SYS IN DD * 


where the underlined are the names of the north, south, and 
plot tapes, in that order, (XR4144 is a copy of XR4026 , 


* See appendix 1. 
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XR4145 of XR4027. If the tape names are changed, change 
them here too.) These are followed by the data cards, 
specifying what plots are to’ be made. 

Each plot is specified by a plot packet, consisting 
of one parameter card, as many coordinate cards as necessary, 
and ending with a delimiter card, in that order. As many 
packets as desired can be put at the end, and the program 
will read one at a time, complete it, and read the next* 

[A time estimate for the job is based primarily on how much 
magnetic tape must be read. If the plot regions are not 
in order in the data cards according to the ordering of the 
tapes (see TAPE FORMAT) north to south, then backspacing will 
occur, but the job will take longer] . 

Contents of the parameter card are specified on that 
list and the card format (Figure 3) . The coordinate cards 
in excess of the number required to specify the plot (either 
1 or 2) are treated as x-ray source positions, and those not 
on the plot region are ignored. (If a star is found within 
EPSLON of the source it is not plotted, but an x-ray symbol 
put next to it in the Table. Each coordinate can be in 
decimal degrees (DD) , radians (RD) , hours, minutes, and 
seconds (HM) , or degrees, minutes, seconds (DM) as specified 
by the letters preceeding it (DD is default) . If it is 
HM, or DM, and negative, the minus sign must go with the 
seconds field only (to prevent trying to read a -0 hours) . 
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The delimiter card must contain asterisks in the first 
two columns. If the next column is blank, there will be no 
title. If it contains an asterisk, a default title is used 
(date and plot number in this job sequence) . If it con- 
tains anything else, col 3-80 are printed as a title. 

The program will compute the center of a limit .specified 
plot on the basis that the geometric center of the specified 
region in the' plot (which will, in general, look like / \ ) 

will be the projection center, so a true picture is obtained. 

The form of the printer table is the same as that in 
the SAO Star Catalog, with the exception that the codes are 
printed as numbers according to the SAO tape code, which is 
appended. (Code 52 is printed explicitly after the index 
number on the left. VV stands for variable in the visual, 

D for double, DV double and variable, etc.) 

The orientation of the plot is always such that vertically 
upwards from the center is NORTH . Right ascension can increase 

i 

from fight to left through the center, or vice versa, depen- 

* 

ding on the INVERT parameter. 

SAOPLOT uses the REREAD routine for data input, and also 
makes use of the PLOT routines available at LNS . The routine 
HACK, which returns the date in a form suitable for printing 
is used, but it can easily be removed. The use of these 
routines has unfortunately made SAOPLOT installation dependent 
and it should not be run on other 360* s (including the Compu- 
tation Center, which does not presently have REREAD) without 



first checking the availability and names of the CALCOMP 
routines and REREAD feature. 

The coordinate card input data have ah original input 
format of DD, since this is the form used in the program. 
Once a specification other than DD has been found on a card, 
this new specification becomes the default (i.e., is used 
until new key letters are found) . This default setting is 
done independently for right ascension and declination; 
changing the default for one has no effect on the other. 

The setting of defaults holds between plots, but is 
always DD for -the first plot of a run. The defaults can 
be easily changed by altering 2 cards in the FORTRAN deck of 
the program as explained in comment cards in the deck. 

The following example may clarify the operation of the 
default system: 

Parameter Card 

1) 5.6634 

2) HM 6 12 -3.8 

3) 7 14 -5.1 

ie-kic 

Parameter Card 

4) 18 23 6.85 

5) DM 220 48 6.231 

*** 


DM ■ 34 26 31.7134 

35 20 14.86 

RD 1.6394 

DD 14.861 
14.256 
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1) declination is in DD by default, right ascension is DM. 

2) declination in HM, right ascension in DM since default 
changed 

3) declination in HM by default, right ascension in RD. 

4) declination in HM since default continues over plot packet 
boundaries, right ascension in DD. 

5) declination in DM, right ascension DD by default. 

Note that negative values under HM and DM must have 
the minus sign on the seconds value. 

111 . The Program INTRP* 

INTRP is a’ Fortran-callable subroutine designed to 
translate the packed 7094 format of the original SAO tapes 
to a format compatible with the 360. Data is passed to it 
by reading a record of the 7094 tapes into a COMMON area 
labeled INREC. The original tapes can be read using the LNS 
routines in the FBSRI system (see the LNS 360 Users Manual). 

They should be read with TRTCH=C. The record length is 2476 
and recording format is U. 

The program returns the translated data in the order and 
form shown in Table 1 in a common area named OUTPAC. 

The form of one record of the old tapes is 50 stars, each 
of which has 11 7094 words of packed data, which comes to 49 1/2 
bytes in the 360. The form of the data for one star in the 
INTRP WORKAREA is shown in Figure 2. The format of the old 
tapes is completely described in the manual for the SAO 

*This is background information and is not used by SAOPLOT. 
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Table I 



Base above 

OUTPAC (bytes) iype 

Quantity 

Bytes 

(total) 

Delta 

0 

real*8 

6 1950 

400 

Alpha 

400 

real*8 

a 1950 

400 

Deilta 0 

800 

real*8 

6 o 

* 400 

W- 

Alpha 0 

1200 

real*8 

a o 

400 

A 

1600 

real*4 

annual y ' (rad) 

200 

B 

1800 

real*4 

annual 0 , (rad) 

200 

C 

2000 

real*4 

r 

annual y (rad) 

200 

D 

2200 

real*4 

annual a (rad) 

200 

E 

2400 

real*4 

photomag 

200 

F 

2600 

real*4 

epoch of <5 0 in 
years-1850 

200 

G 

2800 

real*4 

0 r (rad) 

0 o 

200 

H 

3000 

real*4 

v is mag 

200 

■%7 

V 

3200 

real *4 

epoch of a in 
years~1850 O 

200 

iJ 

3400 

real*4 

0 (rad) 

a 0 

°1950 <rad) 

200 

K 

— 3600 

real*4 

200 

P 

3800 

integer#4 

Spectype (EBCDIC) 

200 

L 

4000 

integer#2 

codes 10&20 

100 

M 

4100 

integer#2 

codes 30 / 40 , 51&52 

100 

N. 

4200 

integer#2 

DM zone & 
codes 60,70 

100 

0 

4300 

integer#2 

DM no. 

100 

Q 

4400 

integer* 4 

Smithsonian book # 
& Source cat. star 

# 200 


4600 



Output common=OUTPAC 
Input common=INREC 
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catalog binary tapes by K. L. Haramundanus , 

A complete listing of INTRP , as well as listings of 
the programs NTAPE and STAPE which contain examples of the 
use of FBSRI and INTRP in the MAIN programs , are in appen- 
dices 3,4, and 5, respectively. 

IV. Tape Format 

One tape contains the northern hemisphere stars; the 
other, the southern. Each tape contains 3600 records, each 
holding stars 'in the area 10n° to 10 (n+1) 0 declination, 

0 . 9K° to 0 .9 (]$+l) ? right ascension, n and K integers. The 
position of a given record r , l.< r <L3 600, tells its field: 
r=400 (8-n) + X since the tapes are ordered similar to the 
SAO catalog tapes: all the 400 records in a 10° band of 

declination are in order of increasing a, and these blocks 
are in order of declination north to south. Data on the 
tapes are at epoch 1950.0. 

The original tapes for the 7094 contain the catalog 
in a, packed format that is convenient for the 7094 but 
rather hard for the 360 to handle. The particular information 
that we require has thus been converted from the old repre- 
sentation to the standard 360 internal floating point format 
(single or double precision as required) . In addition, codes 
that were packed as close as possible in the original tapes 
have been expanded for easier handling in FORTRAN. The 
converted and streamlined version of the catalog was then 
written on two 360 9- track tapes (one each for northern and 
southern hemispheres) . 
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The programs that perform the translation are called 
NT APE and STAPE (two programs are required because of the 
varying number of records per file on the original tapes) . 
They consist of a FORTRAN main program and 4 FORTRAN sub- 
routines to sort the data, switch files on the old tapes, 
write the new tapes, and print a record map of each new 
tape. 

There is also an assembly language program called 
INTRP , which performs the actual translation operations at 
a bit level. INTRP translates all the data, including those 
parts (such as ‘information on fight ascension and declination 
at epoch of original observation) which are not needed by 
SAOPLOT and are therefore not written on the new tapes. 

Anyone having need of this information should see Table 
1 of the output COMMON area of INTRP. Any information that 
is already on the new tapes should be taken from them, since 
they are in a format compatible with the 360. 

V. Parameter Card Quantities 

The quantities on the parameter card (see Figure 3) 
act as follows: 

KEYLETTERS (col. 1-2) tell how the location of the plot 
region is specified: CN, or center, means the next 
card must be a coordinate card containing a legal 
a and 6 which will be used as the center of the 
plot. LM, or limits, means the next two cards 
must be coordinate cards giving the coordinates 
of the lower left and upper right hand corners 
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of the plot, in that order. 

Any further coordinate cards beyond these required, 
in the same packet, will be interpreted as x-ray 
sources, and the first 30 or less will be plotted. 
(Default CN) . 

SCALE (col. 5-10) scale of the plot, in great circle 
degrees per inch at the center, a real number. 
Mandatory for CN plots , calculated from SIZE 
parameter if blanks on LM plot. 

SIZE (col. 15-20) vertical height of the plot, less than 
or equal to the 10" maximum, a real number. (Default 
10 . 0 ) 

MAXMAG (col. 25-30) a real number. Any star from the 

region will be excluded from the plot if its visual 
magnitude is unlisted or greater in numerical value 
(i.e. less bright) than MAXMAG. (Default 20, i.e. no 
stars excluded on this basis) 

MAXNO (col. 35-40) an integer. If the number of stars 
that pass the magnitude test reaches this number, 
a further star is accepted only if it is brighter 
than the dullest star already listed, and it then 
replaces that star in the list. MAXNO can be 
no larger than 30. (Default 30) [note: if MAXNO is 
-1, the star section is skipped, and just the grid 
is drawn] 

EPOCH (col. 70-75) a real number, equal to the year at 
• which the celestial sphere is to be displayed. 
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This plot will be correct for both equinox and 
epoch (i.e., both relative motion of stars, and 
the movement of the vernal equinox will be con- 
sidered) . (Default=1950.0) 

EPSLON (col. 55-65) a real number, expressing a radius 
in degrees. Any stars in the circle of radius 
EPSLON centered on the x-ray source will be entered 
in the list of MAXNO stars as dead entries (i.e., ar$ 
not plotted, but are listed with a source symbol and 
are not subject to elimination if brighter stars are 
found). (Default 10~ 3 ) . (See page 2). 

INVERT (coi. 79) letter I or blank: presence of I means 
a increases from left to right through the plot 
center on the plot. Blank means right to left. 

GRID (col. 45-50) real number, in degrees. A constant 6 
grid line will be drawn on the plot every multiple 
of GRID in declination that falls in the plot region. 
Lines of constant a are also drawn every multiple 
of some "nice" number that makes the grid boxes 
approximately square in the middle. (Default=l. 0) 

[For example, if center a=3.5°, 6=0.7°, GRID=1.0, 
we would have lines at 6=-l, 0, +1, +2, etc., 
a=2 ,3,4] 

TRACE (col. 80) a "T" in this column gives a record of 
which records on the tape were read for this plot. 

SPEC (col. 78) A letter in this column means only stars 
of that spectral type will be included in the list. 
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Example of Plot Packet 


. /*RDCKET_50UECE_3_ 
/ 17j_0_ . . .. 

.475 • 


£ 8.7 


/ 


l c . 


. i 


0 0 0 o'o 0 0 0 0 C 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 o!o 0 0 0 0 (1 0 0 0 0 0 fl 0 0 0 Q 0 OD 0 0 0 0 0 P 

1 2 3 4)5 6 7 8 9 10 1)12 13 14 15 IS 1? 18 19 20 21 22 23 24 2S2S 27 28 29 30 31 32 33 34 35 39 37 38 39 42 43 44 45 <8 47 48 49 59 Si 52 53 54 35 59 57 5* 

1 1 n;i u 1 1 1 1 n 1 1 ii n n 1 1 1 m n 1 1 1 n ;m i i;i 1 1 n 1 1 1 1 11 1 1 1 1 

2 2 2 2j2 2 Z 2 Z 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 1 1 1 1 1 1 1 2 2 2 1 1 l\ll 1 1 1 1 1 1 1 11 1 1 V' 


3 3 3 3 


33333333333333333333333333333333333333333333333 


This packet will produce a plot with center at 
6=15.5°, a=23.7°, and a scale of ,475°/inch [the scale 
on the Palomar plates]. All stars above magnitude 20 
(the default) are selected, up . to 15, after which only 
the brightest are kept. An X-Ray source is plotted at 
6=17.0°, a=22.7°, and all stars in a .1° box around it 
are kept in the list, irrespective of magnitude. The 
Plot is titled "ROCKET SOURCE 3".' A grid of mesh 1° 
(default) is drawn, and the plot will be 10" (default) 
in height. 
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VI . Definitions of Codes in the Smithsonian Catalog 
(From the Smithsonian Tape Description) 


Code 10: Visual Magnitude 


Code 

Photo- 
Visual visual 


Magnitude Source 

Does not appear in source catalog 

Determined by source catalog 

Determined by source catalog or by 
authority in footnote 

Source cited in source catalog introduction 
Source unspecified 
Taken from BD 

Based on Dur chmus terung magnitudes and 
visual estimates 


Taken from AGK 1 


Taken from Cordoba Zones (Resultados) 


Taken from CGA or Cordoba Zones 


Taken from Harvard Publications 

Taken from Harvard or San Luis Photometry 


15 Taken from Henry Draper 

16 Combined magnitude of component stars 

17 Arithmetic mean of maximum and minimum 

magnitudes of a variable star 

Where no figures were reported for magnitudes of a 

variable star, 0.00 was stored on the tapes (i.e. , always 

check code 52 when using magnitudes). When blank, code-0 , 


and field =zeros. 
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Code 20; 

Star 

Number and Footnotes 



Footnote 
without with 

Star Number 



0 

16 

Source catalog only 



1 

17 

Source catalog and BD 



2 

18 

Source catalog and CD 



3 

19 

Source catalog and CPD 



4 

20 

- Cordoba B (Resultados) 

and 

CD 

5 

21 

Cordoba A (Resultados) 

't 

and 

CD 

6 

22 

AGKI and BD 



7 


GC and BD 



8 

24 

Cordoba B (Resultados) 

and 

CPD 

9 


Cordoba A (Resultados) 

and 

CPD 


When blank, code of DM is 0 or 16, footnote is 0 through 9, 
and field is zeros. Footnotes and star numbers are those, 
appearing in. the source catalogs. 

Code 30 : Photographic Magnitude 
Code Source 

0 Does not appear in source catalog 

1 Determined by source catalog 

4 Taken from magnitudes of the CPD and diameters 

of the Cape Astrographic Catalog 

8 Source cited in source catalog introduction 

9 Columbia Contributions Nos. 30 and 31 (Schilt) 
When blank, code is 0, field is zeros. 
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Code 40: Proper Motion 


Code 


Source 


1 

3 

5 


Determined in source catalog 

Determined by comparison of catalog and 
Greenwich AC 

Determined by comparison of catalog and 
AGK 1 


6 Determined by comparison of catalog and 

Greenwich AC on the basis of the smallest 
difference in positions 

w 

8 Determined by comparison of catalog with 

AGK 1 on the basis of the smallest dif- 
ference in positions 

When blank, code is 0, field is zeros. 


Code 51: Spectral Type 

Code Source 

0 Taken from the Henry Draper or no spectrum 
in source catalog 

1 Taken from the HD with M stars reclassified 
by Miss ' Cannon 

2 Classified by G. G. Cillie 

3 Classified by Goedicke 

4 Classified by D. Hoffleit 

5 Classified by M. V. .Maya.ll 

6 Classified by McCormick Observatory 

7 Classified by Nassau and Seyfert 

When no spectra were recorded, the code is 0 and the 
field is BCD blanks. If the spectrum is composite, +++ is 
stored in the field, and the code is 0. 



-14- 


Code 52 : 
Code 
0 
1 
2 

3 

4 

5 

6 
7 


Miscellaneous 

Meaning 

No additional information 

Double star, see source catalog for source 

Double star in Aitken* s Double Star Catalog 

Double star in Burnham's Double Star Catalog 

Variable star in visual magnitude in source 
catalog 

Variable star in photographic magnitude in 
source catalog 

Variable star in both magnitudes 

Both double and variable, in either visual 


or photographic magnitudes 


When blank, code is 0, no field involved 


Code 6 0 : Accuracy of Visual Magnitude 

Listed on tape as 0 or 1; 0 indicates the magnitude was 
reported in source catalog to O^OO; 1, to 0™0 


Code 7 0 : Accuracy of Photographic Magnitude 

Listed on tape as 0 or 1; 0 indicates the magnitude was 
reported in source catalog to 0^0 0 ; 1, to 0™0. 
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VII. Book Numbers for Source Catalogs 

No. Abbreviated Title 


01 

AGK 2 

, vol. 1 



02 

AGK 2 

, vol. 2 



03 

AGK 2 

, vol . 5 



04 

AGK 2 

, vol, 6 



05 

AGK 2 

, vol . 7 



06 

AGK 2 

, vol. 8 



20 

Yale 

Transactions, vol. 

11 


21 

Yale 

Transactions, vol. 

12 

I 

22 

Yale 

Transactions, vol. 

12 

II 

23 

Yale 

Transactions, vol. 

13 

I 

24 

Yale 

Transactions, vol. 

13 

II 

25 

Yale 

Transactions , vol . 

14 


26 

Yale 

Transactions, vol. 

16 


27 

Yale 

Transactions, vol. 

17 


28 

Yale 

Transactions, vol. 

18 


29 

Yale 

Transactions, vol. 

19 


30 

Yale 

Transactions, vol. 

20 


31 

Yale 

Transactions, vol. 

21 


32 

Yale 

Transactions, vol. 

22 

I 

33 

Yale 

Transactions, vol. 

22 

II 

34 

Yale 

Transactions, vol. 

24 


35 

Yale 

Transactions, vol. 

25 


36 

Yale 

Transactions, vol. 

26 

I 

37 

Yale 

Transactions, vol. 

26 

II 

38 

Yale 

Transactions, vol. 

27 


40 

Cape 

Annals, vol. 17 



41 

Cape 

Annals, vol. 18 



42 

Cape 

Annals, vol. 19 




48 


Cape Zone 
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No. Abbreviated Title 

60 Melbourne 3 

61 Melbourne 4 

70 General Catalog 

71 FK3 

74 FK4 





SAMPLE SR0PL0T OUTPUT 


STAR 

NG, alpha 

DELTA 

MAG 

SPECTYPE 

i 

192, 191680 

83,684496 

5,3 

R2 

2 

57,818816 

86 , 488608 

5.8 

F5 

3 

192,-16 1440 

83 i 689552 

5.8 

no 

4 

183; 687200 

87,977056 

6.3 

FG 

5. 

136:318288. 

84 , 386240 

6,3 

FO 

6 

205; 919520 

83 . 003552 

6;. 2 

05 

7 

267-, 300960 

86, 992176 

5.9 

03 

3 

155,674992 

84 . 507856 

5,6 

03 

9 

121, 307200 

84 , 208496 

6,4 

00 

10 

156,291200 

82,814416 

5,3 

F2 

11 

180,540960 

85. 864048 

6.4 

F5 

12 

183,852000 

86,713616 

6, 3 

F2 

13 

27, 203232 

89,028736 

2, 1 

F8 V 

□ 14 

109,460800 

87. 126176 

5, 3 

MO 

15 

267,076000 

86,609552 

4,4 

00 


Figure 

lb. Plotted 

Star Data 




SOURCE 

1 . 

2 . 

3 ’ 


ALPHA 

1 89 , 9998 
276, 7997 
109,4598 


DELTA 

86 . 4999 
87 . 2999 
87, 1259 


Figure lc . Plotted Source Data. 



FORM 




PLOT NUMBER 1 INPUT DATA 



CENTER OF PLO'l — 


DECLINATION 

8 7. 00000 OC 


RIGHT ASCENSION 

180.0000000 


EPOCH 

1950.0 


SCALE 

1 . 000 


SIZE 

lo.o 


MINI HUN MAGNITUDE 

20.0 


MAXIMUM NO. OF STARS 15 

X 

a 

NO. X-RAY -SOURCES 

3 


GR 1 0 • 1.00 


EPS I LON 


5.00E-03 
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* denotes default plot label; blank denotes no plot label; otherwise, the 78 characters 
after the initial double asterisk appear on the printed and plotted output as a label. 



